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An  analytical,  one-dimensional,  steady  state  model  is  employed  to  solve  for  overpotentials  at  the  catalyst 
layers  along  with  the  liquid  water  and  methanol  distributions  at  the  anode,  and  oxygen  transport  at  the 
cathode.  An  iterative  method  is  utilized  to  calculate  the  cell  temperature  at  each  cell  current  density.  A 
comprehensive  exergy  analysis  considering  all  possible  species  inside  the  cell  during  normal  operation  is 
presented.  The  contributions  of  different  types  of  irreversibilities  including  overpotentials  at  the  anode 
and  cathode,  methanol  crossover,  contact  resistance,  and  proton  conductivity  of  the  membrane  are  inves¬ 
tigated.  Of  all  losses,  overpotentials  in  conjunction  with  the  methanol  crossover  are  considered  as  the 
major  exergy  destruction  sources  inside  the  cell  during  the  normal  operation.  While  the  exergy  losses 
due  to  electrochemical  reactions  are  more  significant  at  higher  current  densities,  exergy  destruction  by 
methanol  crossover  at  the  cathode  plays  more  important  role  at  lower  currents.  It  is  also  found  that  the 
first-law  efficiency  of  a  passive  direct  methanol  fuel  cell  increases  as  the  methanol  solution  in  the  tank 
increases  in  concentration  from  1  M  to  3  M.  However,  this  is  not  the  case  with  the  second-law  efficiency 
which  is  always  decreasing  as  the  concentration  of  the  methanol  solution  in  the  tank  increases. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

A  fuel  cell,  in  general,  is  an  electrochemical  energy  device  at 
which  a  fuel  and  oxidant  -  normally  readily  available  oxygen  from 
the  ambient  -  react  at  two  separate  electrodes  to  convert  the 
chemical  energy  into  electrical  energy.  Unlike  conventional  power 
devices,  i.e.  gas  turbines,  steam  turbines,  and  internal  combustion 
engines,  which  are  based  on  certain  thermal  cycles  and  use  irre¬ 
versible  combustion  to  produce  power,  fuel  cells  have  the  potential 
to  provide  more  power  from  a  given  supply  of  fuel  and  oxidant  [1  ]. 
Not  operating  as  a  thermodynamic  power  cycle,  the  notion  of  max¬ 
imum  thermal  efficiency  imposed  by  the  totally  reversible  Carnot 
machine  is  no  longer  valid  for  fuel  cells. 

Among  the  many  different  fuel  cell  types,  direct  methanol  fuel 
cell  (DMFC)  is  categorized  as  a  low-temperature  fuel  cell  with  a 
proton  conductive  polymer  electrolyte  membrane.  DMFCs  are  of 
great  interest  to  provide  power  to  many  portable,  low-temperature 
applications  such  as  laptops  and  cell  phones.  They  offer  advantages 
of  simple  structure  and  instantaneous  charging  time  in  compar¬ 
ison  to  lithium-ion  batteries  [2,3].  However,  there  are  still  some 
barriers  hindering  mass  product  and  commercialization  of  DMFCs, 
including  water  and  methanol  crossover  through  the  membrane, 
and  slow  methanol  oxidation  kinetics  at  the  anode.  Achieving  a 
higher  energy  density  in  DMFCs  is  of  the  highest  priority.  A  higher 
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efficiency,  on  the  other  hand,  results  in  better  and  more  efficient 
utilization  of  a  certain  amount  of  fuel  supply.  Among  all  kinds  of 
analyses,  an  exergy  study  of  DMFCs  offers  a  more  qualitative  knowl¬ 
edge  about  the  exergy  flow  in  and  out,  as  well  as  exergy  destruction 
within  the  system  by  different  kinds  of  irreversibilities. 

Despite  the  large  number  of  exergetic  analyses  conducted  on 
high  temperature  solid-oxide  fuel  cells  for  stationary  power  gen¬ 
eration  [4-7],  there  are  only  few  reported  studies  done  on  the 
exergetic  analysis  of  low-temperature  proton  exchange  membrane 
(PEM)  fuel  cells  [8-12].  Hussain  et  al.  [9]  investigated  the  effect  of 
operating  conditions,  i.e.  temperature,  pressure  and  air  stoichiom¬ 
etry,  on  the  energy  and  exergy  efficiency  of  a  power  cycle  using 
a  PEM  fuel  cell.  They  found  that  the  largest  irreversibility  occurs 
at  the  fuel  cell  stack.  Ay  et  al.  [10]  studied  the  effect  of  operating 
conditions  as  well  as  membrane  thickness  on  the  exergetic  perfor¬ 
mance  of  a  PEM  fuel  cell.  They  reported  a  decrease  in  the  exergy 
efficiency  of  the  cell  with  an  increase  in  the  membrane  thickness.  A 
comparison  was  made  between  DMFC  and  methanol  reforming  of 
PEM  fuel  cell  systems  by  Ishihara  et  al.  [11]  using  exergy  analysis, 
while  the  operating  temperatures  of  both  were  kept  at  80  °C.  They 
noted  that  even  though  DMFCs  do  not  need  a  fuel  treatment  system, 
their  exergy  efficiency  could  not  go  higher  than  that  of  reforming 
PEM  fuel  cells  owing  to  higher  irreversibility  of  methanol  oxidation 
at  the  anode.  Li  et  al.  [12]  investigated  the  exergetic  efficiency  of 
DMFCs,  neglecting  the  effect  of  the  cell  temperature  change.  While 
they  use  an  analytical  solution  to  calculate  the  methanol  crossover 
rate  through  the  membrane,  the  cell  voltage  and  electrochemistry 
are  expressed  by  the  other  experimental  correlations. 
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Nomenclature 

a 

specific  area  for  oxidation/reduction,  m-1 

ACL 

anode  catalyst  layer 

AGDL 

anode  gas  diffusion  layer 

c 

molar  concentration,  mol  nr3 

Cp 

specific  heat  capacity,  J  mol-1  K-1 

CCL 

cathode  catalyst  layer 

CGDL 

cathode  gas  diffusion  layer 

D 

diffusivity,  m2  s-1 

F 

Faraday  constant,  96,478  C  mol-1 

Ex 

exergy,  J 

e 

specific  exergy,  J  mol-1 

g 

gravity,  ms-2 

g 

specific  Gibbs  exergy,  J  mol-1 

Gr 

Grashof  number 

h 

heat  transfer  coefficient,  W  m- 1  K-1 

h 

specific  enthalpy,  J  mol-1 

hfg 

heat  of  vaporization,  J  kg-1 

hm 

mass  transfer  coefficient,  nr2  s-1 

/H+(x) 

protonic  current  density,  Am-2 

JH+  (0)  or  /H+  external  circuit  current  density,  A  nrr2 

'P 

parasitic  current,  Am-2 

rref 

o,02 

reduction  exchange  current  density,  Am-2 

rref 

io,MeOH 

oxidation  exchange  current  density,  A  m-2 

I< 

permeability  of  porous  media,  m2 

k 

thermal  conductivity,  W  nrr 1  K-1 

L 

length,  m 

h" 

molar  flux,  mol  nr2  s-1 

n 

number  of  moles 

Nu 

Nusselt  number 

Hd 

electro-osmotic  drag  coefficient 

p 

pressure,  Pa 

Pr 

Prandtl  number 

q"f 

volumetric  heat  source,  J  nrr3 

R 

reaction  rate,  A  nr3 

Ru 

universal  gas  constant,  J  (mol  I<)-1 

^contact 

Ohmic  contact  resistance,  £2  m2 

5 

entropy,  J  mol-1 1<-1 

S 

liquid  saturation 

Sc 

Schmidt  number 

Sh 

Sherwood  number 

T 

temperature,  K 

U 

thermodynamic  equilibrium  voltage,  V 

X 

coordinate,  m,  or  mole  fraction 

Greek 

a 

transport  coefficient 

y 

reaction  order  of  oxidation/reduction 

£ 

porosity  of  porous  medium 

P 

over  potential,  V 

0 

contact  angle,0 

\l 

viscosity,  kg  nr1  s-1 

P 

density,  kg  nr 3 

o 

interfacial  tension,  N  m- 1 

&  m 

proton  conductivity  of  the  membrane  phase, 
nr1 

<P 

carbon/membrane  phase  potential,  V 

relative  humidity 

Superscripts 

ACL 

anode  catalyst  layer 

amb. 

ambient 

ch 

chemical 

e 

environment 

lim 

limiting  value 

ref 

reference  value 

sat. 

saturated  value 

th-me 

thermo-mechanical 

oo 

values  in  ambient  air 

Subscripts 

a 

anode 

ACL 

anode  catalyst  layer 

AGDL 

anode  gas  diffusion  layer 

amb. 

ambient 

c 

cathode,  or  carbon  phase 

C02 

carbon  dioxide 

cross 

crossover 

CCL 

cathode  catalyst  layer 

CGDL 

cathode  gas  diffusion  layer 

e- 

electron 

g 

gas  phase 

H+ 

proton 

H20 

water 

1 

liquid  phase 

m 

membrane 

MeOH 

methanol 

ox 

oxidation 

o2 

oxygen 

rev. 

reversible 

sat 

saturation  value 

T 

temperature 

V 

vapor 

0 

ambient  condition 

To  the  best  of  authors’  knowledge,  none  of  the  aforemen¬ 
tioned  research  performed  an  exergetic  analysis  for  a  passive  DMFC 
considering  the  temperature  rise  of  the  cell  during  the  normal  oper¬ 
ation.  While  the  operating  temperature  of  an  active  DMFC  is  usually 
kept  fixed  at  a  constant  value,  the  temperature  in  a  passive  DMFC 
is  mainly  governed  by  the  electrochemical  reactions  at  the  cata¬ 
lyst  layers  and  varies  with  the  cell  current  density.  Thus,  providing 
an  exergy  analysis  taking  the  cell  temperature  rise  into  account 
is  inevitable.  Furthermore,  the  analogy  between  heat  and  mass 
transport  offers  the  opportunity  to  have  a  more  comprehensive 
estimation  of  species  concentrations  at  the  cell  boundaries  which, 
in  turn,  results  in  a  more  comprehensive  exergy  analysis  of  the  sys¬ 
tem  in  comparison  to  the  previous  studies  [11,12].  The  objectives 
of  this  study  are  two-fold:  first,  the  fundamentals  of  the  electro¬ 
chemical  kinetics,  transport  phenomena,  and  water  and  methanol 
crossovers  are  developed  through  an  analytical  model;  then,  the 
exergy  analysis  is  utilized  to  provide  a  better  understanding  of  all 
sources  of  entropy  destruction  inside  the  cell. 

2.  Mathematical  models 

The  DMFC,  shown  in  Fig.  1(a),  operates  in  a  completely  pas¬ 
sive  mode  and  is  considered  as  the  thermodynamic  system.  Dilute 
methanol  solution  is  supplied  at  a  constant  concentration  in  the 
fuel  tank.  Methanol  and  water  diffuse  through  the  anode  gas  diffu¬ 
sion  layer  (AGDL)  to  reach  the  anode  catalyst  layer  (ACL)  where  the 
methanol  will  be  oxidized  in  the  presence  of  Pt-Ru  catalysts  to  pro¬ 
duce  carbon  dioxide,  electrons,  and  protons.  The  anode  electrode 
half  reaction  is: 

CH2OH  +  H20  -*  C02  +6H+  +  6e-  (1 ) 
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CGDL 


Fig.  1.  (a)  A  schematic  representation  of  a  passive  DMFC,  and  (b)  protonic  and  electronic  current  density  distributions  at  the  CCL. 


On  the  other  side  of  the  membrane  at  the  cathode,  oxygen 
reaches  the  cathode  gas  diffusion  layer  (CGDL)  by  natural  convec¬ 
tion  from  the  ambient  air.  Then,  it  diffuses  through  the  CGDL  to  the 
cathode  catalyst  layer  (CCL).  At  the  CCL,  oxygen  combines  with  pro¬ 
tons  and  electrons  traveling  through  the  membrane  and  external 
circuit,  respectively,  to  be  reduced  to  liquid  water,  by  the  following 
reaction: 

^02  +  6H+  +6e-  3H20  (2) 

In  order  to  get  the  analytical  solutions  for  the  electrochemistry 
and  distribution  of  the  feeding  species  in  the  system,  i.e.  methanol, 
water,  and  oxygen,  in  conjunction  with  the  cell  temperature,  the 
following  assumptions  are  made: 

•  While  the  electronic  and  protonic  current  densities  change  along 
the  catalyst  layers,  oxygen  and  methanol  concentrations  are 
assumed  to  be  constant  in  these  layers. 

•  The  electrochemical  reaction  rates,  both  at  the  cathode  and 
anode,  are  given  by  Tafel  expressions. 


•  The  electron  conductivity  of  the  carbon  phase  is  assumed  to 
be  much  higher  than  the  proton  conductivity  of  the  membrane 
phase.  As  a  result,  the  carbon  phase  potential,  (pc,  is  assumed  to 
be  constant  at  the  catalyst  layers. 

•  The  convective  transport  of  species  is  neglected  through  the 
entire  domain. 

•  A  constant  concentration  of  methanol  solution  is  supplied  in  the 
fuel  tank. 

•  The  membrane  is  fully  saturated  with  dissolved  water. 

•  Constant,  average  liquid  saturations  are  assumed  for  both  the 
anode  and  cathode. 

•  Temperature  variation  along  the  porous  layers  of  the  cell  is 
neglected,  thereby  providing  a  lumped  method  for  the  energy 
equation. 

2.1.  Electrochemistry 

To  perform  an  exergetic  balance  on  the  DMFC,  providing  an 

analytical  expression  for  the  overpotential  losses  at  the  catalyst 

layers  is  necessary.  In  this  study,  a  similar  procedure  as  that  of 
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[13]  is  taken  to  solve  the  electrochemistry  inside  the  catalyst  lay¬ 
ers. 

The  reduction  reaction  rate  of  oxygen  at  the  CCL  is  estimated  by 
the  following  Tafel  expression  [14]: 


^red  —  (1  ^)^red^oeo2 


0) 


where 


Protons  flowing  through  the  membrane  encounter  internal 
resistance  causing  a  potential  drop  in  membrane  phase  potential 

(0m): 


^resis,H+,mem  —  0m, a  -  0m, c  -  1] 


h+; 


(12) 


where  8m  is  the  membrane  thickness  and  crm  denotes  the  pro¬ 
ton  conductivity  of  the  membrane.  Through  simple  mathematical 
manipulation,  Eq.  (11)  can  be  written  as  follows: 


77c  =  0c  -  0m  -  Uo2  (4) 

Note  that  both  the  electronic  ( /e - )  and  protonic  (/H+ )  current 
densities  can  be  employed  to  calculate  the  overpotentials  at  the 
catalyst  layers.  However,  since  in  this  study  an  infinite  electron 
conductivity  for  the  carbon  phase  is  assumed  (crc  ->  oo),  all  calcula¬ 
tions  are  done  based  on  the  protonic  current  density.  Considering 
the  coordinate  system  shown  in  Fig.  1(b),  the  following  relations 
can  be  written  at  the  CCL: 


dl 


H+ 

dx 


—  ^red 


7j_[+  —  —  crm 


90m 

dx 


0  <  X  <  Lee L 


(5) 


Substituting  Eqs.  (3)  and  (4)  into  Eq.  (5),  one  obtains  an  explicit 
relation  for  the  cathode  overpotential  in  terms  of  cell  current  den¬ 
sity  and  oxygen  concentration  at  the  CCL.  A  detailed  mathematical 
manipulation  is  provided  in  Appendix  A.  The  final  solution  for  the 
cathode  over-potential  is  as  follows: 


i)c  =  0(fH+) '  lnUH+)  -  ln(^o,c)  (6) 

with  the  notation  that  the  bar  sign  refers  to  the  non-dimensional 
parameters  defined  as  follows: 

-  1  -  lccl  •  (1  -  S)  •  ared/£f  2  /  Co2  \  ,/C 

(2amRuT/acF)’  °’c  (2  amRuT/acF)  J  ’ 

V<V)  =  1  +  — ' (7) 
1  +iH+ 


Vcell  =  Uo2  ~  ^MeOH  -  [0c, a  -  0m, a  -  UmcOh] 

-  [Uq2  -  0c, c  +  0m, c]  -  ^resis,H+,mem  (13) 


Taking  into  account  the  contact  resistance  of  the  current  collec¬ 
tor  in  the  cell,  and  following  the  definitions  of  the  cathode  and  the 
anode  overpotentials  given  in  Eqs.  (4)  and  (9),  Eq.  (13)  takes  the 
following  form: 


Vcell  = 


f^MeOH )  -lla~  \rjc\  ~  Jh+ 


(  ^contact  + 


(14) 


where  the  reversible  cell  voltage  is  calculated  by  the  following 
expression  [15]: 

(Uo2-UMeOH)  =  Vrev.  =  Vrev.,298K-(Tcell  ~  298)  X  1.4  X  10“4  (15) 


2.2.  Species  distribution 


It  follows  from  Eqs.  (3)  and  (8)  that  the  oxygen  and  methanol 
concentrations  need  to  be  calculated  at  the  catalyst  layers.  Diffu¬ 
sion  is  deemed  to  be  the  dominant  mechanism  of  species  transport 
in  the  porous  layers  of  the  cell.  Also,  it  is  assumed  that  methanol 
crosses  the  membrane  only  by  diffusion  and  electro-osmosis.  Upon 
using  these  assumptions,  the  methanol  concentration  at  the  ACL  is 
obtained  by  solving  the  following  equation  at  the  anode  side: 


n  dCMeOH 
DMeOH - “ 


fH+  . 

Qp  '  “cross, MeOH 


(16) 


which  can  be  rewritten  as: 


'h-  denotes  the  protonic  current  density  at  the  interface  between 
the  membrane  and  CCL  (at  x  =  0).  Indeed,  this  is  what  is  measured 
in  the  external  circuit  as  the  cell  current  density. 

Likewise,  Tafel  kinetics  is  used  to  describe  the  methanol  oxida¬ 
tion  at  the  ACL,  as  follows: 

Rox  =  «oxfoeMeOH  eXP  (a°x^a  ^ 

where 

Pa  =  0c  —  0m  —  U]\/[eOH  (9) 


Coo  _  f'ACL 
^MeOH  ^MeOH 

f'AGDL 


1 

f^MeOH 


i  rACL 

JH+  ,  n  ^MeOH 

— —  +DMe0H,m— ; - 

or  Lm 


_Ln  ^H+ 

+  nd~jT 


rACL 

'"MeOH 


rACL 

'"MeOH 


+  CACL 
+  lh2o 


(17) 


For  the  sake  of  convenience,  the  following  non-dimensional 
parameters  are  defined: 


rlim 

JMeOH,a 


CMeOHDMeOH  ,lim  CH2QDMeOH 

La gdl  ’  H2°’a  iACDL 


The  final  expression  for  the  anode  overpotential  is  similar  to  that 
shown  in  Eq.  (6): 


p= 


f^MeOH,] 


I'AGDL  t  ^H+ 

-  SMeOH  =  - 


D\ 


MeOH 


rlim 

JMeOH,a 


» £h20  = 


JH+ 

rlim  ’ 
n20,a 


^  =  0(/H+)-ln(/H+)-ln(^o,a)  (10) 

Note  that  for  the  methanol  oxidation,  the  origin  of  the  coordi¬ 
nate,  x  =  0,  starts  at  the  interface  between  the  membrane  and  ACL 
and  ends  at  the  interface  between  the  ACL  and  AGDL.  The  non- 
dimensional  parameters  used  for  methanol  oxidation  are  defined 
in  the  same  fashion  as  those  utilized  for  the  oxygen  reduction  at 
the  cathode  side.  The  cell  voltage,  by  definition,  is  the  difference 
between  the  carbon  phase  potential  of  the  anode  (0c,a)  and  the 
cathode  (0c,c): 


0  = 


rACl 

H20 


Jl 


"MeOH 


H20,ACL 


=  F 


^H20^Me0H 

f'AGDL 


Upon  using  the  above  non-dimensional  parameters,  one  can 
solve  Eq.  (17)  for  the  positive  CA^H,  as  follows: 

rACL 

LMeOH 


(18) 


Vcen  —  0c,  c  —  0c,  a 


(11) 
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with  the  notation  that: 


x  =  1  -  (1  +  6nd)£MeOH  -  (1  +  P)0 


It  should  be  pointed  out  that  6  in  Eq.  ( 1 8 )  is  still  unknown,  which 
requires  solving  for  the  water  concentration  at  the  ACL  Considering 
electro-osmotic  drag  to  be  the  only  mechanism  of  water  crossover 
through  the  membrane,  the  concentration  of  liquid  water  at  the 
ACL  is  obtained  as  the  following: 


/-oo  _  r'ACL 

^h2o  ^h2o 

^AGDL 


dh2o 


7H+  ,  n  7H+ 

CC  d  r 


rACL 

H20 


6  F 


p  f'ACL  .  r'ACL 
'"MeOH  ^H20 


(19) 


Upon  using  the  non-dimensional  parameter,  Z  =  C^q/C^0,  the 
relation  for  Z  appears  as: 


Since  the  local  distribution  of  the  temperature  inside  the  cell 
is  neglected,  by  adding  up  Eqs.  (23)  and  (24),  one  concludes  that 
the  entropy  values  for  electron  and  proton  will  drop  out  from  the 
heat  source  term.  Therefore,  only  the  entropy  change  of  the  overall 
cell  reaction  is  required.  However,  it  worth  it  to  note  that  a  con¬ 
ventional  entropy  system  may  be  used  to  calculate  the  entropy 
of  charged  particles  in  half-cell  reactions  [18-20].  Lampinen  and 
Fomino  [20]  proposed  a  semi-absolute  entropy  system  by  which 
problems  involving  non-charged  and  charged  species  can  be  con¬ 
sidered. 

Generated  heat  at  the  catalyst  layers  is  transferred  to  the  envi¬ 
ronment  by  natural  convection  in  conjunction  with  the  evaporation 
of  condensable  species  in  the  system.  Thus,  the  mass  transport 
has  to  be  taken  into  account.  It  is  a  common  practice  in  thermal 


1  -Z  =  §h2o 


1  + 


12  •  nd  •  (1  +  /3)Z 


(2(1  +P)Z)+  (x+ Vx2  +  (4(l-?MeOH)(l+/8)0) 


(20) 


Since  there  is  no  closed-form  solution  forZ,  Eq.  (20)  can  itera¬ 
tively  be  solved  for  Z.  Note  that  if  the  concentration  of  methanol 
at  the  ACL  was  neglected  compared  to  that  of  water  at  the  ACL,  Z 
would  have  the  following  explicit  expression: 

Z  =  1  -  §h2o(1  +  6nd)  (21) 

Eq.  (21)  can  now  be  used  as  the  initial  guess  to  solve  Eq.  (20). 
Finally,  knowing  the  methanol  crossover  through  the  membrane, 
the  oxygen  concentration  at  the  CCL  is  calculated  as  follows: 

r° °  -  cCCL  J  D 

n  °2  °2  _lH+  ,  3^„  mi 

°2  LCGDL  ~  4F  ^  2  *cross’MeOH  ^ 

2.3.  Cell  temperature 


Providing  a  model  that  can  predict  the  cell  temperature  during 
operation  is  vital  for  an  exergy  study  of  a  passive  DMFC.  However, 
since  the  total  thickness  of  the  five  main  porous  layers  of  the  DMFC 
is  small,  the  spatial  distribution  of  the  temperature  can  be  neglected 
[16,17].  It  is  further  assumed  that  heat  is  only  generated  at  the 
CCL  and  ACL  due  to  electrochemical  reactions.  The  volumetric  heat 
generation  at  the  ACL  due  to  methanol  oxidation  is: 


ASa/6F) 

^ACL  -  (H+  +  T— 


(23) 


where  rja  is  the  anode  overpotential  corresponding  to  the  value 
of  IH+  as  the  cell  current.  At  the  CCL,  oxygen  has  to  be  reduced 
in  the  presence  of  both  electrons  traveling  through  the  external 
circuit  and  electrons  released  by  the  chemical  reaction  of  migrated 
methanol  through  the  membrane.  Considering  this  notion,  the  heat 
source  term  at  the  CCL  is  as  follows: 


^CCL  ~~  (^H+  +  h) 


rjc  -  (T  •  A5C/6F)  T  •  ASt 


f'CCL 


-I? 


6 F  ■  LCc l 


(24) 


where  /p  denotes  the  crossover  current  and  equals: 


Ip  6 F  •  ucross  mcoh 


(25) 


The  Peltier  effect,  (FAS/6F),  for  the  anode  and  cathode  half  reac¬ 
tions  are  calculated  as  the  following: 


/  FAS 


a 


T 
6 F 


(Sco2  +  6SH+  +  6Se- 


SMeOH  -  Sh2o) 


(  FAS 

\~6 T 


T 
6 F 


35h2o  +  ^  Sq2  -  6Sh+  -  6Se- 


(26) 

(27) 


science  to  use  the  analogy  between  heat  and  mass  transport  to  pre¬ 
dict  the  mass  transport  based  on  the  knowledge  of  heat  transport 
phenomenon  as  a  zero-order  approximation.  The  heat  and  mass 
transfer  coefficients  at  interface  (I)  in  Fig.  1(a)  are  taken  from  the 
natural  convection  correlations  on  a  horizontal  surface  facing  up 
[21]. 


Nu  =  0.54(GrPr)°'25,  Nu  = 

Sh  =  0.54(Gr  •  Sc)025 ,  Sh  =  ^ 
Cr_  gPlAp|I3 

ii2 


(28) 

(29) 

(30) 


with  the  notation  that  L  is  the  characteristic  length  and  is  chosen 
as  the  distance  between  the  two  vertical  solid  walls  in  Fig.  1(a). 
The  heat  and  mass  transfer  coefficients  at  the  air  breathing  surface 
(interface  (II))  are  calculated  in  a  similar  fashion  for  a  horizontal 
surface  facing  down.  The  coefficients  for  this  interface  are  as  follows 
[21]: 


Nu  =  0.27(Gr  •  Pr)0'25,  Nu  =  y  (31) 

Sh  =  0.27(Gr  •  Sc)025 ,  Sh  =  ^  (32) 

Evaporation  is  the  other  mechanism  that  cools  the  cell  down. 
Since  the  entire  crossover  methanol  is  consumed  at  the  CCL, 
methanol  evaporation  can  only  occur  at  the  anode  side.  Assum¬ 
ing  that  the  methanol  and  water  are  saturated  at  interface  (I),  the 
molar  flux  of  evaporated  methanol  and  water  at  this  surface  are 
calculated  as: 


^MeOh,v,a  -  ^m,  MeOH  l  LMeOH 


(/ -sat .  _  /~amb.  ^ 
^MeOH  ^MeOH ' 


-  frm,H2o( CSq  ^ H^O ' ) 


^H20,v,a  “  nm,H2' 


(33) 

(34) 


where  superscript  amb.  pertains  to  the  ambient  condition.  While 
the  concentration  of  methanol  vapor  at  the  ambient  is  nearly  zero, 
the  water  vapor  concentration  is  calculated  based  on  the  relative 
humidity  of  the  ambient,  ^amb-,  as  follows: 


T^amb  sat  (jamb.) 

(-amb.  _  _ h2° _ 

H20  -  jamb. 


(35) 


However,  assuming  saturation  condition  for  the  water  vapor  at 
the  cathode  end  of  the  cell  may  result  in  over-predicting  the  evap¬ 
orated  water  from  the  cell  for  small  current  densities.  There  are 
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Anode 


Cathode 


Fig.  2.  Thermal  resistance  network  for  a  passive  DMFC. 


two  sources  of  liquid  water  at  the  cathode  side;  migrated  water 
through  the  membrane  by  the  electro-osmotic  drag,  and  the  water 
generation  due  to  the  reduction  of  oxygen  at  the  CCL.  Both  of  these 
water  sources  are  directly  proportional  to  the  cell  current  density. 
In  this  study,  based  on  numerical  results  [14,16,17],  it  is  assumed 
that  70%  of  the  existing  water  at  the  cathode  is  evaporated  to  the 
ambient,  and  rest  of  that  is  evacuated  in  the  liquid  form. 

Carbon  dioxide  bubbles  at  the  anode  also  carry  some  of  the  evap¬ 
orated  water  and  methanol.  To  have  a  first-order  estimation,  the 
gas  pressure  at  the  anode  has  to  be  calculated.  Numerical  studies 
[14,16,17]  revealed  that  the  liquid  pressure  at  the  anode  as  well 
as  the  gas  pressure  at  the  cathode  is  nearly  equal  to  the  atmo¬ 
spheric  pressure.  Thus,  assuming  a  constant  liquid  saturation,  the 
gas  pressure  at  the  anode  is  obtained  as  follows: 

Pg,a=Pl,a+tfCOS(6>)(|)  J(S)  (36) 

The  Leverette  function,  J(s),  is  a  function  of  liquid  saturation, 
wettability,  and  morphological  properties  of  the  porous  material. 
However,  in  this  study,  Udell’s  correlation  [22]  obtained  from  geo¬ 
logical  porous  media  is  adopted  to  evaluate  J(s): 


J(s)  = 


1.42(1  -s)  —  2.12(1  -s)2  +  1.26(1  -s)3  ifd  <  90° 
1.42s- 2.12s2  +  1.26s3  if0  >  90° 


(37) 


Assuming  that  the  water  and  methanol  vapor  inside  the  carbon 
dioxide  bubbles  are  in  a  saturated  state,  the  following  relation  is 
used  to  calculate  the  partial  pressure  of  C02  inside  the  bubbles: 


Pg,C02  -Pg,a  PMeOH,v  Ph^O.v  (38) 

Knowing  the  rate  of  C02  generation  at  the  ACL,  and  upon 
using  the  ideal  gas  relation,  the  amount  of  evaporated  water  and 
methanol  carried  by  the  C02  bubbles  are  obtained  as  the  following: 


^MeOH ,  v,  bubble,  a 


Jh+  ^MeoH,v 

6P  Pg,C02 


^H20,v, bubble, a 


IH+  Ph2o,v 
6P  Pg,co2 


(39) 

(40) 


The  total  heat  removal  by  the  evaporation  of  water  and 
methanol  in  the  domain  are  as  follows: 


^H20,evap.  (^H20,v,bubbl,a  ^H20,v,a  4"  ^H20,v,c)  ^fg,H20  (41) 

^MeOH.evap.  =  (^MeOH,v,bubbl,a  +  ^MeOH,v,a^fg,MeOH  (42) 

Fig.  2  displays  the  thermal  resistance  network  involving  heat 
generation  and  heat  removal  at  both  the  anode  and  cathode.  Note 
that  an  iterative  manner,  with  an  appropriate  under- relaxation,  has 


to  be  taken  to  get  the  converged  value  of  the  cell  temperature  at 
each  individual  current  density  of  the  cell. 

2.4.  Exergy  balance 

The  second  law  of  thermodynamics  puts  restrictions  on  the  effi¬ 
ciency  of  transforming  one  form  of  energy  to  another.  As  a  result, 
exergy  is  a  property  that  quantifies  the  potential  for  maximum 
work.  Unlike  energy  which  is  conserved  in  every  process,  exergy 
can  be  destroyed  by  the  irreversibilities  or  be  transferred  to  and 
from  thermodynamic  systems.  Depicted  in  Fig.  3(a)  is  a  schematic 
of  all  the  components  with  their  corresponding  pressure  and  tem¬ 
perature,  in  a  DMFC  during  the  passive  operation.  As  shown,  there 
are  a  total  of  three  components  coming  into  the  system  and  six  dif- 


the  environment  at  T0 


Fig.  3.  A  schematic  representation  of  (a)  all  components  coming  into  the  system 
(DMFC)  and  going  out,  with  their  corresponding  pressure  and  temperature,  and  (b) 
the  reactor  used  to  introduce  the  chemical  exergy  of  liquid  methanol 
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Table  1 

Exergy  reference  environment. 


Thermo-mechanical  exergy  reference  environment 

T0  298  K 

Po  101,325  Pa 

Chemical  exergy  reference  environment 

Gas  phase  species 

Molar  fraction 

n2 

0.7560 

02 

0.2034 

H20  (vapor) 

/( 1/7),  i/r  =  70% 

co2 

0.0003 

ferent  species  going  out.  While  the  cell  produces  some  electrical 
work  and  exchanges  heat  to  the  environment,  some  of  the  input 
exergy  is  destroyed  due  to  irreversibilities.  Looking  at  Fig.  3,  the 
exergy  balance  of  a  DMFC  may  be  written  as  follows: 


EEx  mass,  in  —  EEx  1 


-mass, out  + 


FxWork  +  ^  ^  FXheat  +  FXi0St 


(43) 


with  the  notation  that: 

y^£xmass,in  =  if1"  *  ^)MeOH,l  +  *  ^)o2  +  '  ^)h20,1  (44) 

^  ^£Xmass,out  =  (h  •  £)h20,1  +  '  ^)c02  +  '  ^)h20,v 

+  (”"  •  e)MeoH,v  (45) 


(4S| 

For  a  reacting  system,  the  total  exergy  consists  of  both  thermo- 
mechanical  and  chemical  exergies.  Calculating  the  chemical  exergy 
requires  quantifying  the  molar  fractions  of  all  available  species  in 
the  exergy  reference  environment,  which  is  presented  in  Table  1. 
Assuming  that  the  partial  pressure  of  the  species  of  interest  in  the 
reference  environment  is  P0,  the  thermo-mechanical  exergy,  by 
definition,  is  the  maximum  theoretical  work  obtainable  from  that 
species  of  property  of  (T,  P)  as  it  comes  into  thermal  and  mechanical 
equilibrium  with  the  reference  environment.  The  general  equation 
for  thermo-mechanical  exergy  is  expressed  as  the  following: 

=  &  ~  T°~skP  ~  Cho  ~  ToSoWpo)  (47) 

where  the  ideal  gas  relations  are  used  to  calculate  the  enthalpy  and 
entropy  difference  for  all  gaseous  species  in  the  domain. 

On  the  other  hand,  chemical  exergy  is  the  maximum  theoret¬ 
ical  work  that  can  be  produced  by  a  virtual  Carnot  machine  in 
which  a  substance  of  interest  enters  at  (T0,  Po)  and  completely 
expands  or  reacts  to  produce  the  environmental  components  at 
(T0,  XjP0),  where  x2-  denotes  the  molar  fraction  of  that  component 
in  the  environment.  It  should  be  pointed  out  that  for  the  liq¬ 
uid  state  of  the  condensable  species  in  the  system,  i.e.  water  and 
methanol,  the  chemical  exergy  is  defined  as  the  maximum  theoret¬ 
ical  work  that  a  Carnot  machine  can  produce  as  these  components 
enter  the  machine  at  (T0,  Po)  and  exit  at  (T0,  xzP0),  where  x2-  is  the 
molar  fraction  of  the  vapor  state  of  corresponding  species  in  the 
environment. 

2.4 A.  Exergy  flow  into  the  system 

Since  methanol  enters  the  cell  at  (T0,  Po)  in  the  liquid  form,  the 
contribution  of  thermo-mechanical  exergy  is  zero.  Flowever,  the 
chemical  exergy  is  not  zero.  Fig.  3(b)  shows  how  one  mole  of  liq¬ 
uid  methanol  reacts  with  the  environmental  component  of  oxygen 
to  produce  other  corresponding  environmental  components.  The 


maximum  work  obtainable  from  this  system,  chemical  exergy,  is 
equal  to  the  Gibbs  energy  difference  of  the  reactants  and  products, 
as  follows: 


eMeOH,in  =  (eCh)(T0,P0)  =  C&R  “  ^Sp) 


&MeOH,l  +  2Sc°2  ~SC02  ~  2gH20,\ 


-I  (T0,Po) 


(48) 


+  PuTq  In 


*0, 


^AH20,v  aC02  y 


where  the  Gibbs  energy  of  oxygen  at  the  atmospheric  conditions  is 
zero. 

Since  the  water  in  the  exergy  reference  environment  exists  in 
a  vapor  form,  liquid  water  carries  some  chemical  exergy  into  the 
system,  as: 


eH20,l,in 


_  (gth-me)^^)  +  (ech )(T0,P0) 

=  0  +  f  [&h2o,i  -  #H2o,v](T0)p0)  +  Ku T0  In 


1 

ve 

aH20,v 


(49) 


Note  that  since  the  oxygen  enters  the  system  as  an  environmen¬ 
tal  component,  both  its  thermo-mechanical  and  chemical  exergies 
are  zero. 


2.4.2.  Exergy  flow  out  of  the  system 

Exergy  of  carbon  dioxide  going  out  of  the  system  at  the  anode 
is  obtained  as  follows: 


eco2,a,out  =  (e1 1  )(r0,xco,,aPo)  +  (e  W^o) 


(Cp, co2(TCeii  -  To))  -  To  (Cp, co2  In  (^ 


To  ) 


(50) 


-Ru  In 


fPc O^c 
V  Po 


RUT0  In 


Exergy  of  C02  at  the  cathode  is  calculated  in  a  similar  fashion. 
Water  vapor  also  takes  some  exergy  out  of  the  system.  The  exergy 
of  water  vapor  at  the  anode  is  obtained  as: 


^H20,v,a,out  — 


-Ru  In 


/PFI2O,  v,  a 
V  Po 


RuT0  In 


H20,v  /  J 


(51) 


Exergy  taken  out  of  the  system  by  methanol  vapor  at  the  anode 
can  be  calculated  in  a  same  fashion  as  of  Eq.  (48)  with  the  exception 
of  using  the  vapor  phase  Gibbs  energy  for  the  methanol.  The  chem¬ 
ical  exergy  of  liquid  water  exiting  the  cathode  side  at  (Tcell,  Plc) 
is  the  same  as  that  of  liquid  water  at  the  anode,  and  therefore,  Eq. 
(49)  can  be  used  for  this  purpose.  Flowever,  the  thermo-mechanical 
exergy  has  to  be  calculated  using  Eq.  (47)  assuming  an  incompress¬ 
ible  liquid.  Assuming  a  constant  liquid  saturation  at  the  cathode,  the 
cathode  liquid  pressure  is  calculated. 

3.  Results  and  discussion 

Physicochemical  properties,  as  well  as  cell  geometries  and  oper¬ 
ating  conditions,  are  given  in  Tables  2  and  3,  respectively.  All 
electrochemical  and  physical  parameters  in  the  present  model 
are  calibrated  using  the  polarization  experimental  data  of  [23]. 
As  shown  in  Fig.  4(a),  a  good  agreement  between  the  experi¬ 
mental  data  and  results  obtained  by  the  present  analytical  model 
is  seen  for  two  different  methanol  solutions  in  the  tank.  Since 
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a  non-isothermal  model  is  presented,  a  validation  for  the  pre¬ 
dicted  temperature  is  necessary.  Depicted  in  Fig.  4(b)  is  the 
comparison  between  the  cell  temperature  reported  by  [23]  and 
that  obtained  by  the  present  model.  It  should  be  pointed  out 
that  since  [23]  used  a  tank  with  a  finite  capacity,  the  cell 
temperature  never  reaches  the  steady  state  value.  The  experi¬ 
mental  cell  temperature  data  comprises  the  cell  temperature  after 
the  first  15  min  of  each  cell  load  being  applied.  Temperatures 
obtained  by  the  present  model  are  within  2%  of  the  experimental 
data. 

Of  particular  interest  in  the  DMFC  is  the  crossover  of  water 
and  methanol  through  the  membrane.  There  are  two  mechanisms 
of  methanol  crossover  considered  in  this  study;  electro-osmosis 
and  diffusion,  as  shown  in  Fig.  5(a).  While  the  diffusive  mecha¬ 
nism  of  methanol  crossover  decreases  monotonically  versus  cell 
current  density,  the  electro-osmotic  drag  of  methanol  shows  an 
increase  at  low  currents  and  then,  after  reaching  a  certain  pick 


point,  starts  decreasing.  The  slope  of  the  increase  in  the  electro- 
osmotic  drag  of  methanol  at  low  currents  increases  as  the  methanol 
solution  in  the  tank  becomes  more  concentrated.  Summation  of 
these  two  mechanisms  is  the  total  methanol  crossover  flux  through 
the  membrane  which  is  shown  in  Fig.  5(b).  While  for  1  M  and 
2M  solutions  in  the  tank  the  total  flux  of  methanol  monotoni¬ 
cally  reduces  as  the  current  increases,  methanol  crossover  for  3  M 
solution  first  shows  an  increase  at  low  current  densities  and  then 
decreases  after  reaching  a  maximum  value.  The  underlying  ratio¬ 
nale  for  this  behavior  is  that  the  methanol  concentration  at  the 
ACL  decreases  as  the  cell  current  increases  (Eq.  (17)).  Although  the 
electro-osmotic  crossover  of  methanol  is  linearly  proportional  to 
the  cell  current  density,  the  reduction  of  CMeoH,ACL  is  the  dom¬ 
inant  factor  in  decreasing  the  total  methanol  crossover  as  the 
cell  current  density  increases  for  methanol  feeding  concentrations 
lower  than  3  M.  The  variation  of  the  methanol  crossover  versus 
cell  current  density  for  different  methanol  solutions  are  in  good 


Table  2 

Physicochemical  properties. 


Parameters 

Symbols 

Value 

Unit 

Ref. 

MeOH  in  water  I 
Diffusivities,  J 

9 MeOH 

1.58  x  10-9e0,02623(T~298)(£. s)T 

m2  s'1 

[26] 

liquid  MeOH  in  Nafion  I 

<■  DMe0H,m 

4.9  x  10-10e[2436('333_^](£,i5)'r 

m2  s'1 

[26] 

Diffusivities,  gas 

db,o2 

T  1.823 

1775  *  10-5  (273.15)  Mi-S))T 

m2  s'1 

[27] 

Electron  conductivity  in  carbon  phase 

Proton  conductivity  in  membrane  phase 

°c 

Gm 

4000 

3 

fl  _1m_1 

[28] 

r 

Lp,l,H20 

/— 2.667  x  1(T7  T3  +  2.68  x  KT4  T2  +\ 

\  -8.92  x  10-2  T  +  1.4023  x  101  )  *  18 

J  mol'1  K'1 

Cp,l,MeOH 

79.5 

J  mol'1  K'1 

Cp,H2/Ru 

/  -1.812  x  10“12  74  +  5.521  x  10“9  T3  \ 

\-5.8  x  10“6  T 2  +  2.67  x  10“3  T  +  3.057 ) 

- 

Heat  capacity  { 

Cp,o2  /Ru 

(  2.156  x  10-12T4- 6.764  x  10"9r3+  \ 
\7.055  x  10“6  T2  -  1.878  x  10“3  T  +  3.626/ 

- 

[21] 

fitted 

Cp,co2/Ru 

/2.002  x  10“9  T3  -  6.607  x  10~6  T2\ 

\  +8.735  xl0"3r  + 2.401  J 

- 

Cp,g,H2o/Ru 

/  0.807  x  10-12  T4  -  2.964  x  10 ~9T3+  \ 
\4.152  X  10“6  T2  -  1.108  X  10“3  T  +  4.07 ) 

- 

Cp,g,MeOH 

f— 2.0833  x  10~7  T 3  +  1.895  x  10"4  T2\ 

\  -5.385  xl0-2T  +  6.187767  J  * 32 

J  mol1  K1 

Absolute  entropy 
(1  atm ,  298  K) 


Gibbs  function  of 
formation 
(1  atm,  298  K) 


Liquid  phase 


Gas  phase  < 


Liquid  phase  ■< 


Gas  phase  < 


Electro-osmotic  drag  coefficients 
Henry  law  constant  for  methanol 


c° 

s  l, MeOH 
c° 

bl,H20 

S°C02 


S02 


SH2 


dl.MeOH 


9l,H20 


9co2 


9v,h2o 

9 1, MeOH 

nd 

ku  MpHH 


126.8 

69.95 

213.685 

205.033 

128.7 

-166290 

-237180 

-394380 

-228590 

-162140 

2.5 

0.096e[ao4511(r_273)] 


J  mol 1  K1 
J  mot1  K'1 
J  mol1  K'1 
J  mol'1  K'1 
Jmof1  K1 
J  mol'1 

J  mol'1 

J  mol'1 
J  mol'1 
J  mol'1 

atm 


[1] 


[27] 

[27]  and 
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Table  2  ( Continued ) 


The  saturation  pressure  of  pure  water 

pure 

ioavrw 

Latent  heat  of  evaporation 

hfg,MeOH 

(20  <  T(  °C)  <  80 ) 

hfg,H2o 

Cell  thermodynamic  voltage 

Vrev.,2%  K 

Transfer  coefficient  of  anode 

CCa 

Transfer  coefficient  of  cathode 

ac 

Anode  exchange  current  density 

jref 

lo,MeOH 

Cathode  exchange  current  density 

}ref 

lo,02 

| 

f  aox 

Specific  area  < 

1  ared 

Anode  reference  concentration 

rref 

Ll,MeOH 

Cathode  reference  concentration 

rref 

lo2 

Surface  tension  ^ 

%0 

(10  <  T(  °C)  <  75 ) 

aMeOH 

Liquid  water  density' 

Pi,h2o 

Liquid  methanol  density 

PlMeOH 

-2,1794  +  0.02953(7  -  273.15) 


-9.1837  x  10-s(r  -  273.15)2 
+1.4454  x  10-7(r  -  273.15)3 

atm 

[26] 

(  -1,292 r3  +  76.25T2  -  \  1Q-i 

V14558.3T  + 1.2182  x  107/ 

Jkg1 

[26] 

/ -1,25  x  KT1T3  +  4.999T2  A  1Q-i 

1 23,6  xl03r  +  2,5009  xlO7/ 

Jig1 

1.216 

V 

[15] 

1. 

- 

Assumed 

1.1 

- 

Assumed 

[35570/  1  1\1 

94.25el  Ru  ^53"rJ] 

Am2 

[26] 

[73200/  1  IX] 

0.04222el  Ru  '^"rj] 

A  m2 

[26] 

106.1 

m 

5000 

m1 

Assumed 

10000 

mol.m 3 

Assumed 

renv. 

lo2 

mol.m 3 

Assumed 

/9.3  X  103r3  -  3.96T2\  x  10-7 

V  -1363.87  +  756325  / 

Nml 

[19] 

/3.3xl0~2r2  -774.5nyin-7 
l  +240041.67  rW 

Nm1 

Fitted 

/6.9094  -  2.0146  x  10_5(7  -  273)\ 
-5.9868  x  10"6(r  -  273)2  + 
exp  2.5921  x  10_8(7  -  273)3  - 
9.3244  xl0-n(7-  273)4  + 

\  1.2103  x  10-13(7- 273)5  / 

Kg.m3 

[28] 

244.4  X  0.224(  ^  513^) 

Kgm3 

[29] 

agreement  with  that  has  been  reported  by  Meyers  and  Newman 
[24]. 

In  DMFCs,  methanol  crossover  is  strongly  interconnected  with 
the  water  crossover.  DMFCs  using  highly  concentrated  methanol 
can  achieve  a  low  methanol  crossover,  if  the  water  crossover  is 
limited  [14].  Since  in  this  study  only  the  electro-osmotic  drag  of 
water  through  the  membrane  is  considered,  an  increase  in  water 
crossover  is  observed  as  the  cell  current  increases.  If  the  molar  frac¬ 
tion  of  liquid  water  at  the  ACL  is  approximated  to  be  unity,  as  is  the 
case  for  most  analytical  solutions  such  as  [12,13],  there  would  be 
no  difference  in  water  crossover  for  different  methanol  concen¬ 
trations  in  the  tank.  Flowever,  Fig.  6(a)  illustrates  that  the  molar 
fraction  of  liquid  water  at  the  ACL  is  less  than  unity,  especially  for 
lower  cell  current  densities  and  higher  methanol  concentrations  in 


Table  3 

Cell  geometric  dimensions  and  operating  parameters. 


Parameters 

Symbols 

Value 

Unit 

Anode  gas  diffusion  layer  thickness 

6agdl 

0.26 

mm 

Anode  catalyst  layer  thickness 

<$ACL 

0.02 

mm 

Membrane  thickness  (Nation  115) 

6mem 

0.125 

mm 

Cathode  catalyst  layer  thickness 

6ccl 

0.02 

mm 

Cathode  gas  diffusion  layer  thickness 

6cgdl 

0.26 

mm 

Permeability 

K 

2  x  10-12 

m2 

Porosity 

s 

0.6 

- 

Contact  angle 

0 

120 

- 

Anode  liquid  saturation 

Sa 

0.8 

- 

Cathode  liquid  saturation 

Sc 

0.1 

- 

the  tank.  Fig.  6(b)  depicts  that  the  water  concentration  at  the  ACL 
reduces  as  the  current  density  increases  due  to  higher  electrochem¬ 
ical  consumption  and  water  crossover  through  the  membrane.  As 
expected,  the  water  concentration  at  the  ACL  is  lower,  when  the 
methanol  solution  in  the  tank  is  more  concentrated. 

As  previously  noted,  there  are  three  components  carrying  chem¬ 
ical  exergy  into  the  system;  oxygen  at  the  cathode,  and  liquid 
water  and  methanol  at  the  anode.  Flowever,  since  it  is  assumed 
that  the  oxygen  at  the  cathode  is  supplied  by  the  environment,  the 
total  exergy  of  oxygen  always  equals  zero.  Fig.  7  shows  the  total 
exergy  flow  of  liquid  water  and  methanol  at  the  anode,  which  is 
completely  contributed  by  the  chemical  exergy.  As  expected,  the 
chemical  exergy  of  liquid  methanol  flowing  into  the  system  is  sev¬ 
eral  orders  of  magnitude  greater  than  that  of  liquid  water.  It  is 
observed  that  the  total  flow  of  exergy  increases  as  the  cell  current 
density  increases,  which  is  a  direct  result  of  a  higher  mass  flow  rate 
and  cell  temperature  at  larger  current  densities.  Exergy  supplied  by 
the  liquid  methanol  to  the  system  is  one  order  of  magnitude  larger 
than  that  by  the  liquid  water. 

On  the  other  hand,  there  are  a  total  of  six  components  flow¬ 
ing  out  of  the  system,  among  which  methanol  vapor  at  the  anode 
carries  the  highest  exergy  out  of  the  cell.  Methanol  evaporation, 
either  from  the  top  of  the  anode  tank  (interface  (I))  or  into  the  C02 
bubbles,  strongly  depends  on  the  cell  temperature  and  methanol 
concentration.  As  shown  in  Fig.  8,  the  exergy  of  methanol  vapor  for 
3  M  methanol  solution  is  larger  than  that  for  1  M.  While  the  C02  pro¬ 
duction  and  its  exergy  flow  at  the  cathode  are  proportional  to  the 
methanol  crossover  through  the  membrane  and  tend  to  decrease  at 
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Fig.  4.  (a)  Polarization  curve,  and  (b)  the  cell  temperature  versus  current  density, 
for  two  different  methanol  solutions  in  the  tank. 


high  current  densities,  exergy  flow  of  C02  at  the  anode  monotoni- 
cally  increases  with  current  density,  mostly  owing  to  the  higher  rate 
of  production  and  cell  temperature.  The  flow  of  exergy  for  vapor 
and  liquid  water  are  also  shown  in  Fig.  8. 

Fig.  9  shows  the  exergy  balance  of  the  cell,  as  described  in  Eq. 
(43),  for  two  different  methanol  solutions  in  the  anode  tank.  At 
open  circuit  voltage  of  the  cell  (7  =  0  Am-2),  there  is  some  exergy 
flow  in  due  to  the  liquid  water  and  methanol  crossovers  through 
the  membrane  and,  also,  methanol  evaporation  from  the  anode. 
Flowever,  the  entire  migrated  methanol  through  the  membrane  is 
chemically  reacted  with  water  and  oxygen  at  the  CCL  without  pro¬ 
ducing  any  useful  power.  Heat  exergy  loss  to  the  environment  is 
directly  related  to  the  cell  temperature  and,  as  seen,  increases  as 
the  cell  current  goes  up.  Though  its  value  is  relatively  small  com¬ 
pared  to  the  exergy  flow  in  and  overall  exergy  loss,  the  heat  exergy 
loss  to  the  environment  is  considerable  in  comparison  with  the 
network  exergy. 

The  cell  performance  is  degraded  due  to  exergy  destruction— a 
measure  of  the  magnitudes  of  the  irreversibilities  present  during 
operation.  There  are  multiple  sources  of  irreversibilities  that  exist 
during  the  normal  operation  of  a  DMFC,  including;  phase  change  of 
condensable  species,  friction,  mixing  of  multiple  components,  elec¬ 
tric  resistance,  and  the  electrochemical  reaction.  The  total  exergy 
loss  for  a  DMFC  during  operation  versus  current  density  for  three 
different  methanol  concentrations  in  the  tank  is  shown  in  Fig.  9. 
However,  more  detailed  information  about  the  lost  exergy  and  the 
contribution  by  each  of  the  sources  are  provided  in  Fig.  10,  where 
the  term  “Other”  refers  to  the  exergy  loss  due  to  the  phase  change, 


(a) 

la. 

*3 


a 


x  IQ-3  X  IQ'3 


X  IQ'3 


Fig.  5.  (a)  The  electro-osmotic  drag  and  diffusive  transport  of  methanol  through 
the  membrane,  and  (b)  total  methanol  crossover  flux  for  three  distinct  methanol 
concentrations  in  the  tank. 


mixing  of  the  components  in  the  system,  and  friction.  Exergy  loss 
due  to  the  electrical  resistance  in  the  membrane  and  the  contact 
resistance  is  relatively  small,  and  defined  as  follows: 


£*lost,res.  —  ^H+2  f  ^contact  +  — — 
\  om 


(52) 


Half  reactions  at  the  anode  and  cathode  are  in  equilibrium  condi¬ 
tions  when  no  external  overpotential  voltage  exists.  This  condition 
results  in  a  zero  net  flux  of  reactants  and  products  in  the  reac¬ 
tion  of  interest.  However,  during  the  normal  operation  of  a  DMFC, 
external  overpotentials  are  required  at  both  the  anode  and  cath¬ 
ode  to  force  the  half  reactions  to  proceed  in  a  desired  direction  at 
an  acceptable  rate.  These  cause  the  electrochemical  irreversibilities 
inside  the  cell,  either  at  the  anode  or  the  cathode,  and  are  defined 
as  follows: 


£*lost,a,elec.chem.  —  ^h+  '  (53) 

£*lost,c,elec.chem.  —  ^h+  '  rlc 

The  values  of  the  anode  and  cathode  overpotential  exergy  losses 
are  directly  related  to  the  corresponding  transfer  coefficients,  ared 
and  c^ox,  which  are  calibrated  by  the  experimental  data  of  [23].  As 
shown  in  Fig.  10,  the  exergy  losses  due  to  the  anode  and  cathode 
overpotentials  are  very  significant  at  higher  cell  current  densities. 
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Fig.  6.  Variation  of  (a)  molar  fraction  of  water  at  the  ACL,  and  (b)  liquid  water 
concentration  at  the  ACL,  versus  cell  current  density. 

However,  variation  of  exergy  loss  contributed  by  the  methanol 
crossover  is  different  than  that  by  the  overpotentials.  While  the 
methanol  crossover  exergy  loss  monotonically  decreases  as  the  cell 
current  density  increases  for  the  case  of  1  M  methanol  concentra¬ 
tions  at  the  tank,  an  initial  increase  followed  by  a  final  reduction  is 
observed  when  3  M  methanol  concentration  is  fed  at  the  anode. 


Fig.  7.  Exergy  flow  of  liquid  water  and  methanol  into  the  system. 


3B 

£ 


V(A.m-2) 
(a)  1M 


V(A.m-2) 
(b)  3M 


Fig.  8.  Exergy  flow  out  of  the  system  by  six  different  components  versus  cell  current 
densities,  when  the  methanol  concentrations  in  the  tank  are  (a)  1  M,  and  (b)  3  M. 


One  of  the  objectives  of  the  present  study  is  to  use  the  energy 
and  exergy  concepts  to  assess  the  effectiveness  of  the  DMFC  during 
normal  operation.  To  this  end,  the  first-law  and  second-law  effi¬ 
ciencies  are  introduced.  The  first-law  efficiency  of  any  conversion 
device  is,  by  definition,  the  ratio  between  useful  energy  outputs 
to  the  total  energy  inputs.  In  the  case  of  a  DMFC,  the  maximum 
first-law  efficiency  is  obtainable  when  all  Gibbs  free  energy  can  be 
converted  to  electrical  energy,  and  is  defined  as: 

AC 

'H  I,  max  —  T77  (54) 


where  the  lower  heating  value  of  fuel  is  used  throughout  this  study. 
The  maximum  first-law  efficiency,  ^/max,  is  independent  of  the  cur¬ 
rent  density  and  takes  the  values  of 0.967  and  0.973  when  liquid  and 
vapor  methanol  is  respectively  supplied  at  the  anode  as  fuel.  How¬ 
ever,  the  real  first-law  efficiency  is  much  lower  than  these  values 
and  is  defined  as  the  ratio  of  useful  work  to  the  energy  input: 


^useful 

(AH)  x  (7/6F) 


(55) 


Though  the  first-law  efficiency  gives  some  insight  about  the 
effectiveness  of  a  DMFC,  it  cannot  assess  the  effectiveness  of  energy 
resource  utilization.  To  this  end,  the  work-second-law  efficiency  is 
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Fig.  9.  Exergy  balance  for  the  cell  at  different  cell  current  densities  when  the 
methanol  concentrations  in  the  fuel  tank  are  (a)  1  M,  and  (b)  3  M. 


introduced  as  follows: 


^H.work  -  W°^ £  (56) 

/^CAmass,in  —  /  y£*mass,out 

In  the  same  fashion,  the  heat-second-law  efficiency  -  a  mea¬ 
sure  of  assessing  exergy  loss  by  heat  transfer  from  the  cell  to  the 
environment  -  is  defined  as: 


77ll, Heat 


£*Heat 


y^^xmass.in  —  X/  *mass>out 


(57) 


Shown  in  Fig.  11  are  the  first-  and  second-law  efficiencies  ver¬ 
sus  cell  current  density  for  three  different  methanol  concentrations 
in  the  anode  fuel  reservoir.  While  an  increase  in  the  first-law  effi¬ 
ciency  is  observed  for  higher  methanol  solution  concentrations,  the 
work-second-law  efficiency  increases  as  the  methanol  concentra¬ 
tion  in  the  tank  decreases.  It  should  be  pointed  out  that  the  first-law 
efficiency  does  not  show  any  optimum  condition  for  the  opera¬ 
tion  of  a  DMFC,  while  according  to  the  work-second-law  efficiency, 
there  is  always  an  optimum  condition  for  the  cell  at  a  specific  cur¬ 
rent  density.  These  optimum  current  densities  are  1210, 1030,  and 
810  Am-2  for  3, 2,  and  1  M  methanol  solutions,  respectively,  where 
the  corresponding  second-law  efficiencies  are  0.122,  0.162,  and 
0.224.  The  lower  work-second-law  efficiency  corresponding  to  a 
higher  methanol  solution  is  attributed  to  the  higher  irreversibil¬ 


-  Over-potential -Anode 


(a)  1M 


V(  A.m-2) 
(b)  3M 


Fig.  10.  Different  types  of  irreversibilities  present  in  a  DMFC  versus  cell  current 
when  the  methanol  concentrations  in  the  fuel  tank  are  (a)  1  M,  and  (b)  3  M. 


ities  associated  with  the  methanol  crossover  (Fig.  10).  Another 
conclusive  point  in  Fig.  1 1  is  the  relatively  small  heat-second-law 
efficiency  (^n.Heat)  compared  to  the  work-second-law  efficiency 
( *7n,work)>  mainly  due  to  the  passive  operation  of  the  cell  in  which  the 


Fig.  11.  First-  and  second-law  efficiencies  of  a  passive  DMFC  versus  cell  current 
density  for  three  methanol  concentrations  in  the  tank. 
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V  (A.m-2) 


Fig.  12.  The  effect  of  ambient  temperature  on  the  second-law  efficiencies  when  1  M 
methanol  solution  is  fed  in  the  fuel  tank. 

cell  temperature  is  driven  by  the  internal  electrochemical  reactions 
rather  than  an  external  heat  source. 

The  assessment  of  an  energy  conversion  system  using  the 
second-law  analysis  is  highly  dependent  on  the  definition  of  the 
exergy  reference  environment  (dead  state).  Depicted  in  Fig.  1 2  is  the 
effect  of  the  exergy  reference  ambient  temperature  on  the  work- 
and  heat-second-law  efficiencies  when  1  M  methanol  concentra¬ 
tion  is  fed  at  the  anode  fuel  reservoir.  As  seen  in  this  figure,  the 
colder  reference  ambient  temperature  results  in  better  second-law 
efficiency. 


Fig.  Al.  Proton  current  density  distribution  in  the  CCL  for  different  cell  current 
density. 


Fig.  1  (b)  displays  the  solution  domain  and  the  coordinate  system 
of  solving  the  cathode  overpotential.  Upon  combining  Eqs.  (3),  (4) 
and  (5)  and  using  the  non-dimensional  parameters  of  Eq.  (A-l ),  one 
arrives  at  Eq.  (A-2): 

x  -  X  1  -  1  _  rj  _  rj 

Lccl  {^crmRuT/acF)  I*  {RuT/(xcF)  77* 

=  (A-l) 


4.  Conclusions 

The  exergy  analysis  of  a  passive  DMFC  provides  more  detailed 
information  about  different  types  of  irreversibilities  present  inside 
the  cell.  The  following  conclusions  are  made  upon  reviewing  the 
above  results: 

•  Methanol  crossover  in  conjunction  with  the  anode  and  cathode 
overpotentials  contributes  the  largest  portion  of  the  irreversibil¬ 
ities  inside  the  cell.  While  the  methanol  crossover  contribution  is 
mostly  significant  for  lower  current  densities,  the  exergy  loss  due 
to  overpotentials  are  more  important  at  higher  current  densities. 

•  The  work-second-law  efficiency  of  a  passive  DMFC  is  mono- 
tonically  decreases  as  the  methanol  solution  becomes  more 
concentrated.  In  other  words,  a  lower  concentration  of  methanol 
in  the  tank  always  offers  a  better  work-second-law  efficiency. 

•  Since  the  cell  temperature  of  a  passive  DMFC  is  only  affected  by 
the  electrochemical  reactions  at  the  catalyst  layers  rather  than  an 
external  heat  source,  the  heat-second-law  efficiency  of  the  cell  is 
low  compared  to  the  work-second-law  efficiency. 

•  Cell  operating  at  lower  ambient  temperature  offers  a  better 
second-law  efficiency. 
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Appendix  A.  Overpotentials  at  the  cathode  catalyst  layer 

Both  protonic  and  electronic  current  densities  can  be  used  to 
derive  the  overpotentials  at  the  anode  and  cathode.  However,  it 
is  more  convenient  to  use  the  protonic  current  density,  since  the 
electron  conductivity  of  the  carbon  phase  is  chosen  to  be  infinity. 


92/h+  9(/h+)2 

dx2  dx 

Rewriting  Eq.  (A-2),  the  following  is  concluded: 

+  j£+(x)  =  -Rred(x)  +  /2  +  (x)  =  -Rred(0)  +  /2  +  (0) 


(A-2) 


(A-3) 


where  ~IH+  (0)  denotes  the  current  density  at  x  =  0,  which  is  identical 
with  the  measured  current  at  the  external  circuit.  Also,  other  non- 
dimensional  parameters  are: 


Rred(x)  =  R0,c-exp(7)c(x)) 

f'CCL  '  (1  —  s)  '  ared^oeo2 
Rq.c  =  - rr - 1 — 


Tef\^c 


(Co2/C£) 


(A-4) 


The  right  hand  side  of  Eq.  (A-3)  has  to  be  negative  [13].  Consid¬ 
ering  that  the  CGDL  is  impermeable  to  the  proton  current  density, 
the  solution  of  Eq.  (A-3)  is 


/H+(*)  =  yKred(o)-/2+ 


Eq.  (A-5)  shows  an  implicit  relation  between  /H+(0)and  Rred(0). 
For  a  specific  /H+( 0),  infinite  solutions  for  Rred(0)  exist.  However, 
since  there  is  just  one  realistic  solution  for  the  overpotential,  and 
Rred(0)  has  to  be  finite,  the  smallest  Rred(0)  is  the  physical  solution. 
Fig.  Al  shows  the  protonic  current  distribution  at  the  CCL  for  dif¬ 
ferent  cell  current  densities.  As  seen,  an  increase  in  the  cell  current 
density  results  in  a  more  nonlinear  distribution  of  current  density 
at  the  CCL.  A  linear  distribution  assumption  of  current  density  is 
valid  only  for  low  current  densities. 

Unfortunately,  an  explicit  relation  between  /H+(0)  and  Rred(0) 
cannot  be  obtained.  However,  Kulikovsky  [25]  solved  Eq.  (A-5), 
at  x  =  0,  for  small  and  large  values  of  /H+( 0)  using  the  following 
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approximation 


References 


X  •  tan(A)  = 


A2  A  «  1 
A  A  »  1 


(A-6) 


An  explicit  expression  for  cell  current  density  in  terms  of  non- 
dimensional  cathode  overpotential  is  obtained  for  two  extreme 
cases,  as  follows: 

(  R0,c  ■  exp(j)c(0))  7h+(0)  «  1 

/H+(°)={^exp(5f))  /H+(0)»l  ^7) 

Now,  overpotential  can  easily  be  written  in  term  of  current  den¬ 
sity  of  the  cell. 


-  (m=f  ln(/H+(0))  -  ln(R0,c)  /H+(0)  «  1 
mUJ  \  2  •  ln(/H+(0))  -  ln(Ro,c)  /H+(0)  »  1 


(A-8) 


Looking  at  Eq.  (A-8),  one  concludes  that  the  only  difference 
between  the  overpotentials  at  two  extreme  non-dimensional  cur¬ 
rent  densities  is  the  first  term  at  the  right  hand  side  of  Eq.  (A-8). 
In  order  to  have  the  above  relation  for  the  entire  range  of  non- 
dimensional  current  density,  a  fitting  function  is  used  as  the 
following  [25]: 


rjc(O)  =  WH+(0))  •  ln(/H+(0))  -  ln(R0,c)  (A-12) 


where 


VK/h+(0))  =  1  + 


1 

l+/H+(0) 


(A-13) 


The  same  procedure  is  taken  for  the  anode  overpotential.  It 
should  be  pointed  out  that  all  the  non-dimensional  parameters 
have  to  be  redefined  based  on  the  anode  parametric  values.  Also, 
note  that  in  the  text  for  the  sake  of  convenience,  ~IH+  is  used  in 
replace  of/H+(0). 
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